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Kvll.l channels are critical for the maintenance of a normal heart rhythm. The flow of potassium ions through 
these channels is controlled by two voltage-regulated gates, termed "activation" and "inactivation," located at op- 
posite ends of the pore. Crucially in Kvll.l channels, inactivation gating occurs much more rapidly, and over a 
distinct range of voltages, compared with activation gating. Although it is clear that the fourth transmembrane 
segments (S4), within each subunit of the tetrameric channel, are important for controlling the opening and clos- 
ing of the activation gate, their role during inactivation gating is much less clear. Here, we use rate equilibrium free 
energy relationship (REFER) analysis to probe the contribution of the S4 "voltage-sensor" helix during inactivation 
of Kvll.l channels. Contrary to the important role that charged residues play during activation gating, it is the 
hydrophobic residues (Leu529, Leu530, Leu532, and Val535) that are the key molecular determinants of inactiva- 
tion gating. Within the context of an interconnected multi-domain model of Kvll.l inactivation gating, our REFER 
analysis indicates that the S4 helix and the S4-S5 linker undergo a conformational rearrangement shortly after 
that of the S5 helix and S5P linker, but before the S6 helix. Combining REFER analysis with double mutant cycle 
analysis, we provide evidence for a hydrophobic interaction between residues on the S4 and S5 helices. Based on a 
Kvll.l channel homology model, we propose that this hydrophobic interaction forms the basis of an intersubunit 
coupling between the voltage sensor and pore domain that is an important mediator of inactivation gating. 
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INTRODUCTION 

Kvll.l channels play an important role in repolariza- 
tion of the cardiac ventricular action potential and are 
therefore key regulators of a normal electrical rhythm 
in the heart (Sanguinetti and Tristani-Firouzi, 2006) . 
Like other voltage-gated potassium (Kv) channels, Kvll.l 
channels are comprised of four subunits, each with six 
transmembrane segments (S1-S6) (Vandenberg et al., 
2012). On the periphery of each subunit are the S1-S4 
helices, which form the voltage-sensing domain. The 
central conduction axis of the channel is lined by the 
four pore domains, formed by the S5 and S6 helices 
separated by a short pore helix and selectivity filter. Two 
voltage-regulated gates, located at either end of the 
pore domain, control the flow of potassium ions through 
Kvll.l channels. The activation gate, formed by the in- 
tracellular ends of the four S6 segments, is slow to open 
and close. In contrast, the inactivation gate, formed by 
the four selectivity filter segments at the extracellular 
side of the pore, has rapid gating kinetics (Smith et al., 
1996; Spector et al, 1996; Wang et al, 2011). These 
unique kinetics properties make Kvll.l channels ideally 
suited for determining the duration of the plateau 



phase of the cardiac action potential (Vandenberg 
etal., 2012). 

Alterations to the voltage-dependent inactivation gat- 
ing of Kvll.l channels can have profound effects on 
the duration of cardiac repolarization and hence dura- 
tion of the Q_T interval on the body surface electrocar- 
diogram. Mutations that shift the voltage dependence 
of inactivation to more depolarized potentials result in 
reduced inactivation leading to greater current flow 
and shorter repolarization intervals (Brugada et al., 2004; 
Sun et al., 2011), whereas mutations that shift the volt- 
age dependence in the hyperpolarized direction have 
the opposite effect and result in lengthening of repolar- 
ization intervals (Zhao et al., 2009). Both of these scenar- 
ios result in a marked increase in the risk of arrhythmias 
and sudden cardiac arrest. 

Although both the activation and the inactivation 
gates of Kvl 1 . 1 channels are voltage regulated, the two 
gating processes occur over different voltage ranges 
(Wang et al., 1997) and exhibit different sensitivities 
to amino acid mutations (Zou et al., 1998; Sanguinetti 
and Xu, 1999; Piper et al., 2005; Clarke et al., 2006), 
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temperature (Vandenberg et al., 2006), and extracellu- 
lar [Ca 2+ ] (Johnson et al., 1999). Although the first 
three to four positive charges of the S4 helix have been 
clearly established as the voltage sensor for Kvll.l acti- 
vation (Subbiah et al., 2004; Zhang et al., 2004; Piper 
et al., 2005), the role of the S4 helix during Kvll.l 
channel inactivation is relatively unknown. Even though 
the two gating processes are quite distinct, there is no 
reason why the S4 voltage sensor could not undergo two 
distinct movements and so underlie both functions. 
The unusual gating kinetics of Kvll.l channels allows 
for separation of the rapid open-to-inactivated gating 
transition from the much slower closed-to-open transi- 
tion, thereby permitting an examination of the S4 helix 
contribution to inactivation gating, independent of the 
S4 helix motion that transfers positively charged resi- 
dues across the membrane electric field (gating cur- 
rent) and is coupled to activation gate opening. 

In a previous study (Wang et al., 2011), we showed 
that the dynamics of inactivation gating in Kvl 1.1 chan- 
nels was amenable to rate equilibrium free energy re- 
lationship (REFER) analysis. REFER analysis compares 
the effect of systematic perturbations to a protein on 
the free energy changes in kinetic (AAG + ) and thermo- 
dynamic (AAGo) relationships between two stable end 
states to gather information on the "structure" of the 
transition state (Fersht et al., 1992; Fersht, 2004, 2008) 
and, as a consequence, on the relative timing of domain 
motions during the transition between two stable end 
states (Grosman et al., 2000; Cymes et al., 2002; Zhou 
et al., 2005; Auerbach, 2007; Purohit et al., 2007) . When 
applying REFER analysis to proteins, mutations of in- 
dividual residues can be used to gain molecular-level 
detail about the conformational changes taking place 
(Grosman et al., 2000; Fersht, 2004; Zhou et al., 2005; 
Purohit et al., 2007; Wang et al., 2011). The slope of a 
REFER plot, that is the slope of A AG 1 versus AAGo, is 
defined as the <J>-value, and for this reason, in the pro- 
tein-folding literature REFER analysis is usually refer- 
red to as $-value analysis (Fersht et al., 1992; Fersht, 
2004, 2008). ^-Values that are between 0 and 1 are in- 
dicative of the relative time point at which the mutation 
perturbed the interconversion between the two stable 
states, where <J>-values close to 0 indicate that the muta- 
tion has affected late steps in the reaction, and <J>-values 
close to 1 indicate that the mutation perturbed an early 
step in the reaction (Grosman et al., 2000; Cymes et al., 
2002; Zhou et al., 2005; Auerbach, 2007; Purohit et al., 
2007; Wang et al., 2011). In our previous study (Wang 
et al., 2011), we found that mutations to charged resi- 
dues in the S4 region produced relatively modest per- 
turbations to inactivation gating. In this study, we used 
REFER analysis to undertake a more comprehensive 
analysis of the role of the S4 helix during inactivation 
gating of Kvll.l channels. Our results show that hydro- 
phobic residues within the S4 domain are crucial in 



mediating inactivation gating of Kvll.l channels, inde- 
pendendy of the role that the charged residues of the 
S4 helix play during channel activation. 

MATERIALS AND METHODS 



Molecular biology and channel expression 

Female Xenopus laevis frogs were purchased from Nasco, and all 
animal procedures were approved by the Garvan and St. Vincent's 
Animal Ethics Committee (AEC 08/34, 1 1/37) . Frogs were anaes- 
thetized by immersion in 0.17% wt/vol tricaine, and stage V or VI 
oocytes were removed via a small (~7-8-mm) abdominal incision. 
Oocytes were either mechanically separated using forceps, or lobes 
of oocytes were digested with 1 mg/ml collagenase A for 2 h, and 
then stored in ND96 solution (mM: 96 NaCl, 2 KC1, 1.8 CaCl 2 , 
1 MgCl 2 , and 5 HEPES, pH 7.5) supplemented with 10 ug/ml 
gentamicin and 2.5 mM pyruvate. 

Kvl 1 . 1 cDNA (provided by G. Robertson, University of Wisconsin- 
Madison, Madison, WI) was subcloned into a pBluescript vector 
containing the 5' untranslated region (UTR) and 3' UTR of the 
Xenopus (3-globin gene (provided by R. Vandenberg, University of 
Sydney, Sydney, Australia) . Mutations were introduced using the 
Quikchange method (Agilent Technologies) and confirmed by 
DNA sequencing. Plasmid DNA was linearized using BamHI-HF 
(NEB), and the in vitro transcription of cRNA was performed 
using the mMessage mMachine kit (Ambion). cRNA encoding 
WT or mutant (mut) Kvll.l channels was injected into oocytes, 
which were then incubated at 18°C for 18-72 h before electro- 
physiological recording of current. 

Electrophysiology 

Oocytes were impaled with glass capillary micropipettes that had 
access resistances in the range of 0.3 to 1.0 MO. Currents were 
recorded from oocytes using a Geneclamp500B two-electrode 
voltage-clamp amplifier interfaced to a PC via a Digidata 1440 
(Molecular Devices). Signals were filtered at 2 kHz and digi- 
tized at 5-10 kHz. Voltage-clamp protocols are shown as insets 
in Fig. 1 A. In all protocols, the cells were held at —90 mV before 
any voltage steps. After each experiment, the offset potential was 
checked, and if it exceeded ±5 mV, the recordings were discarded. 
Data acquisition and analysis were performed using pCLAMP 
(version 10.2; Molecular Devices), Excel (Microsoft) and Prism 
(version 6; GraphPad) software. All parameter values were calcu- 
lated as mean ± SEM for n experiments, where n denotes the 
number of different oocytes studied for each construct. 

Kinetic analysis 

The first step in performing a REFER analysis is to determine 
the kinetics of both the forward (open-to-inactivated) and reverse 
(inactivated-to-open) gating transitions of WT and mutant chan- 
nels. This is relatively straightforward for Kvll.l channels be- 
cause the rapid inactivation gating transitions can be easily 
separated from the open/closed transitions (Smith et al., 1996; 
Spector et al., 1996; Wang et al., 2011). Rates of inactivation were 
measured using a triple-pulse protocol (see Fig. 1 A, i) in which 
cells were depolarized to +40 mV (or +80 mV for some mutant 
channels) for 500 ms to activate and inactivate channels, and then 
stepped to —90 or —110 mV for 10 ms to allow recovery from 
inactivation, before stepping to voltages in the range of +100 to 
— 60 mV to monitor the onset of current inactivation. A single- 
exponential function was fitted to the decaying portion of the 
corresponding current trace (see Fig. 1 B, i) . Reverse rates of in- 
activation gating (inactivated-to-open, k rcc ) , normally termed "re- 
covery from inactivation," were measured from two-step voltage 
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protocols in which cells were depolarized to +40 mV (or +80 mV) 
for 500 ms to activate and inactivate channels, and then stepped 
to voltages in the range of —40 to —180 mV to allow channels 
to rapidly reopen before slowly closing (Fig. 1 A, ii) . The initial 
rapid component of the "hooked" tail current at each voltage 
(depicted by black lines in Fig. 1 B, ii) represents the inactivated- 
to-open-state transition, whereas the slower component repre- 
sents the open-to-closed-state transition. Therefore, the faster 
time constant obtained from a double-exponential function fit- 
ted to the hooked tail current represents the rate of recovery 
from inactivation. 

For a bimolecular reaction with voltage-dependent rates, such 
as Kvll.l channel inactivation, the observed rate (hi hs ) at any par- 
ticular voltage (V) is equal to the sum of the forward (£j nact ,v) and 
reverse rates (^cv) measured at that voltage: 

Kla,V = hnatf.V + Kec.V (1) 

At 0 mV, &obs,o is therefore a combination of hmaft + &rec,o- However, 
at very positive voltages, k, sc is extremely small, and hence k„ bs » 
kinaa- Conversely, at very negative voltages, Ai nact is extremely small, 
and so « k TCC . An extrapolation of the linear portions of fto bs at 
very positive or very negative voltages (dotted lines in Fig. 1 D) 
can then be used to derive the unidirectional forward (A inact0 ) and 
reverse (Arec.o) rate constants at 0 mV, respectively. A plot of the 
logarithm of fcbs.v results in a characteristic chevron phenotype 
(Fig. 1 D, i), indicating a reaction that is dominated by a single 
transition step (Jackson and Fersht, 1991). The equilibrium con- 
stant for inactivation at 0 mV (-Ke q ,o) was calculated by: 

K,q,0 ~ ^inact.O / (2) 

Perturbations to inactivation gating for each mutant were mea- 
sured as the Alog(Ae q0 ) of the mutant channel compared with WT: 

Alog(^ q , 0 ) = log(^ q , 0 ) mut - log^.oW (3) 

A positive shift in Alog(^ eq> o) is equivalent to a hyperpolarizing 
shift in the voltage dependence of inactivation, whereas a nega- 
tive shift in Alog(-K,. q? o) is equivalent to a depolarizing shift in the 
voltage dependence of inactivation. In addition, the limiting 
slopes of the fits to the chevron plots were used as a direct mea- 
sure of the voltage dependence of the forward and reverse transi- 
tion pathways. 

For a given mutation, the <t>-value is calculated by comparing 
changes, versus WT, in the energetics of the transition state 
(AAG*) relative to changes in the free energy difference between 
the ground states (AAG 0 ) . For a two-state reaction, the transition 
state energy, AG*, is given by: 

AG t =-RT.(ln(/ £iW )+ ln(A)), (4) 

where A is the "pre-exponential factor" (Price and Dwek, 1996). 
The effect of a mutant on the transition-state energy is given by 
the equation: 

AAG* = -RT. (ln(A inactWT ) - \n(k hmriM ) ) , (5) 

i.e., we assume that the mutation does not alter the "pre-exponential 
factor" for the reaction. Similarly, the effect of a mutant on the 
equilibrium free energy is given by: 

AAG 0 = -RT. (ln(^ iWT ) - ln(^ mut ) ) . (6) 

A <t>-value can be calculated from the ratio of the change in the 
logarithm of the unidirectional forward rate constant, log(A inact0 ), 



relative to the change in the logarithm of the equilibrium constant, 
log(Ae q , 0 ): 

O = AAG* / AAG 0 = Alog^o) / Alog(^ >0 ). (7) 

The $-value calculates to what degree a mutation-induced 
perturbation in the equilibrium constant for inactivation, which 
represents the energy difference between the open and inacti- 
vated ground states, is caused by a change in the forward transi- 
tion rate (see Fig. SI). The <I>-value then reflects at what stage during 
the native reaction pathway the mutated residue experienced a 
change in its environment, as reflected by the change to the ener- 
getics of gating. To derive an accurate $-value, a mutation must 
cause a sufficient perturbation to the equilibrium constant, i.e., 
the denominator in Eq. 7. We, both previously (Wang et al., 201 1 ) 
and in Fig. S2 B, and others (Cymes et al., 2002; Fersht and Sato, 
2004) have established a Alog(i^e q0 ) of greater than or equal to 
±0.5 as the cutoff criterion to derive an accurate $-value. Informa- 
tive <l>-values are those that lie between 0 and 1 (Cymes et al., 
2002; Wang et al., 2011), as shown in Fig. SI. There are three 
scenarios that result in an invalid $-value: (1) the mutation af- 
fects the transition state but not the ground states (i.e., a catalytic 
mutation), resulting in an infinite $-value (Fig. SI F); (2) the 
mutation produces opposite effects on the transition state com- 
pared with the equilibrium between the two stable ground states, 
resulting in a negative $-value (Fig. SI G); and (3) the mutation 
produces a greater effect on the transition state than on the 
ground state, resulting in a $-value >1 (Fig. SI H) . One benefit to 
REFER analysis is that these uninformative mutations can be iden- 
tified because they result in $-values that lie outside of the range 
of 0 and 1 . It is, however, also possible that a mutation could affect 
two or more independent processes that, by coincidence, cancel 
each other out to give a <5-value between 0 and 1. Because this 
coincidental combination is unlikely to occur with a diverse range 
of mutations at a given residue, a better estimation of the $-value 
is derived from families of mutations at that position (Grosman 
et al., 2000). The other advantage of using families of mutants is 
that if a single mutant gives a very different result than that of 
other mutants at the same residue, it is immediately very obvious 
(see Fig. S2 A). 

Double mutant cycle analysis 

To test for an energetic coupling between two residues in different 
helices, we combined REFER analysis with double mutant cycle 
analysis. In this approach, individual mutations in the S4 helix, 
which significantly perturb inactivation (measured as Alogi^q o > 
±0.5 log units relative to WT), were combined with perturbing 
mutations in the S5 helix (Ile560, Leu564, or Ile567) or the S5P 
linker (Asp591). If an energetic coupling exists between two 
residues, the perturbations caused by the individual mutations 
(AAG mutl and AAG mut2 ) would not be additive when combined in 
the double mutant (AAG mutl+mut2 < AAG mut i + AAG mut2 ) . Conversely, 
if the perturbations caused by the two single mutants were addi- 
tive when combined in the double mutant (AAG mut i +mut2 = AAG mut i 
+ AAG mut2 ) , the two residues are not energetically coupled. 

Mean (X), SEM, standard deviation (S), and number of obser- 
vations (n) for AAG mutl + AAG mut2 were calculated by summing 
the corresponding values of Alog(Ae q>0 ) from the two individual 
mutant channels (mutl and mut2). An unpaired, unequal vari- 
ance Welch's t test was used to perform statistical analysis be- 
tween the measured perturbation caused by the double mutant 
(AAG mutl+mut2 , termed "X double ") and the summed values of the two 
single mutants (AAG mm i + AAG mut2 and X mm i +mut2 ). For the statisti- 
cal test, t-values were calculated as: 

t = -^(double) - -^(mutl+mut2) ' ^"(^(double) _ -^(mutl+mut2))' (8) 
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where the variance (S(X(d 0 ubie) X(muti+mut2) ) ) was calculated by; 



^ v v ,S(double) 2 ^ S(mut\ + mut2) 2 ^ 

' '"' n(double) n(mut\ + mttil) 

Degrees of freedom (df) was calculated by: 

— (^double / "double) + (^mutl+mut2 / n mutl+mut2 ) / (( (^double /"double) / 
("double - W+HS-ltJ /"muH+muts) 2 / ("mrtl+mutS ~ J ) ) ) • 

(10) 
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Figure 1. Measurement of inactivation kinetics for WT and 
V535S Kyi 1.1 channels. (A) Voltage protocols used to measure 
the rates of onset of inactivation (i) and the rates of recovery from 
inactivation (ii) , with key regions for the measurement of current 
highlighted in black. (B and C) Representative families of current 
traces measuring the rates of onset of inactivation (i) and rates of 
recovery from inactivation (ii) for WT (B) and V535S (C) Kvll.l 
channels. Highlighted in bold are current traces at 0 mV (i) and 
— 130 mV (ii) to aid comparison. (D) Chevron plots of the loga- 
rithm of the observed rates for onset of (open squares) and recov- 
ery from (closed squares) inactivation for WT (i) and V535S (ii) 
channels, plotted against voltage. Solid black lines are a fit of Eq. 1 
(see Materials and methods), whereas the solid gray line in (ii) 
indicates for WT channels to aid comparison. Dashed lines 
indicate the derived unidirectional rate constants for the onset 
of (knaay) and recovery from (hecy) inactivation, with values at 
0 mV indicated by arrows. The equilibrium constant for inactiva- 
tion (Ae q ) at 0 mV was calculated by: = Ajnaco/Areco (Eq- 2 in 
Materials and methods). 



p-Values were calculated using t and df values. P < 0.05 was consid- 
ered significant. 



Voltage dependence of activation 

Isochronal activation curves were measured according to stan- 
dard tail current analysis (Vandenberg et al., 2012). From a hold- 
ing potential of —90 mV (or —120 mV for some mutants), the 
membrane potential was stepped to voltages in the range of — 60 
to +50 mV (the precise range varies from —120 to +100 mV de- 
pending on the mutant studied) for 4 s, before a step to —70 mV 
(or —120 mV) to record tail currents. Tail currents were normal- 
ized to the maximum current amplitude (7„ im ) and fitted with a 
Boltzmann expression: 



(ID 



where ///„,„ is relative current, 1^.5 is the voltage of half-current 
activation, V, is the membrane potential, and k is the slope factor. 
Alternatively, data were fitted with the thermodynamic form of 
the Boltzmann expression: 



(&G 0 -z e EF) 

l + e 



(12) 



where A Go is the work done at 0 mV, z g is the effective number of 
gating charges moving across the membrane electric field (E) , Fis 
Faraday's constant, R is the universal gas constant, and T is the 
absolute temperature. From Eq. 12 we can calculate the effect of 
mutations on changes in chemical potential energy: 



AAG„ = AG„ 



-AG„ 



(13) 



Homology structure generation 

The Kvll.l homology model, generated using Swiss PdbViewer 
(Guex and Peitsch, 1997) and optimized using SWISS-MODEL 
Workspace (Arnold et al., 2006; Bordoli et al., 2009) , was based on 
the crystal structure of a Kvl. 2/2.1 channel chimera (Long et al., 
2007) , according to the sequence alignment shown in Fig. S3. 

Online supplemental material 

Fig. SI shows a series of energy diagrams that explain the theory 
of REFER ($-value) analysis when applied to the study of Kvll.l 
channel inactivation gating. Fig. S2 demonstrates the variability in 
<l>-values derived from mutants that perturb Kvll.l channel inac- 
tivation with varying severity. Fig. S3 is a sequence alignment 
between Kvll.l and Kvl. 2/2.1 channels, which was used to gen- 
erate our Kvll.l channel homology model. Table SI summa- 
rizes the kinetic parameters of inactivation for each S4 mutant 
investigated, and Table S2 summarizes the activation parame- 
ters for each S4 mutant investigated. The online supplemental 
material is available at http://www.jgp.org/cgi/content/full/ 
jgp.201310975/DCl. 



RESULTS 

REFER analysis can be applied to any process where 
there are two ground states separated by a single domi- 
nant transition state, and where the forward and reverse 
rates can be well described by a single-exponential 
process (Auerbach, 2007). C-type inactivation gating in 
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Figure 2. Scanning mutagenesis of 
the Kvll.l channel S4 helix. (A-C) 
Shifts, relative to WT, in log(Ae qi0 ) val- 
ues for alanine (A), serine (B), and 
tryptophan (C) mutagenesis scans of 
S4 residues spanning from Gly522 to 
Lys538. Data are presented as means 
± SEM for n = 3-20 cells (see Table SI) . 
Dashed lines indicate Alog(K eq0 ) > ±0-5 
log units, which has been shown previ- 
ously to be the minimum perturbation 
required to derive an accurate $-value 
(Cymes et al., 2002; Wang et al., 2011). 
Mutations that fulfill this minimum re- 
quirement are indicated by closed bars, 
whereas mutations with Alog(^ qi0 ) < 0.5 
log units are indicated by open bars. 
*, mutant channels that failed to ex- 
press or expressed poorly. 



Kvll.l channels is very rapid, so the forward (open-to- 
inactivated) and reverse (inactivated-to-open) inactiva- 
tion gating transitions can easily be separated from the 
much slower activation (closed-to-open) or deactivation 
(open-to-closed) gating transitions (Smith et al., 1996; 
Spector et al., 1996; Wang et al., 2011). Although it 
is possible that Kvll.l channels may also transition 
into an inactivated state from a preopen closed state 
(Kiehn et al., 1999), this transition is likely negligible 
when compared with the open-to-inactivated transition 
pathway (Bett et al., 2011). Furthermore, the probability 



of a closed-to-inactivated transition can be minimized 
using voltage protocols that ensure that the majority of 
channels are in the open conformation before initiat- 
ing inactivation. 

The specific voltage protocols used to measure rates 
of inactivation (hma, measured at voltages ranging from 
—20 to +80 mV) and rates of recovery from inactivation 
(&rec; measured at voltages from —60 to —160 mV) are 
shown in Fig. 1 A, with the resulting currents for WT 
Kvll.l channels (Fig. 1 B) or V535S channels (Fig. 1 C) 
shown below. For WT, the region of current representing 
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Figure 3. Families of mutations at important 
Kvll.l channel S4 residues. (A and B) Shifts in 
log(K cq fi ) , relative to WT, for families of mutations 
at charged (A) or hydrophobic (B) residues that 
were identified as important determinants for in- 
activation gating from our mutagenesis scans of 
the S4 helix. Mutations are indicated by single- 
amino acid code on the y axis. Data are presented 
as means ± SEM for n = 4-14 cells (see Table SI). 
Mutations that cause a significant perturbation 
to inactivation, measured as Alog(-Ke q! o) > ±0.5 
log units relative to WT, are indicated by closed 
bars, whereas mutations with Alog(-Ke q 0 ) < ±0.5 log 
units are indicated by open bars. *, mutant chan- 
nels that failed to express or expressed poorly. 
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the onset of inactivation (Fig. 1 B, i) or the recovery 
from inactivation (Fig. 1 B, ii) at a given voltage (V) 
were well described by a single-exponential function 
and were used to measure ^ nac t,v and Kecv, respectively. 
A plot of the observed rates (-Kobs! see Eq. 1 in Materials 
and methods) against voltage generates a "chevron" 
phenotype, which is indicative of a transition domi- 
nated by a single step (Fig. 1 D, i). The extremities of 
each arm of the chevron plot are linear as they reflect 
the unidirectional rate constants of the forward or re- 
verse transition at the extreme voltages. To calculate an 
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Figure 4. Mutation of S4 residues affects inactivation at a mid- 
way point during the transition pathway. (A) REFER plots of the 
log forward unidirectional rate constant, log(A inact0 ), against the 
log equilibrium constant, log(A^ Ji o) . f° r inactivation, in families of 
mutations at individual residues: i, Leu529; ii, Leu530; iii, Leu532; 
iv, Val535. (B and C) REFER plots for all mutations at Leu529, 
Leu530, Leu532, and Val535 combined (B) or all S4 residue mu- 
tations (C). Data are presented as means ± SEM for 3-14 cells. 
The slope of the linear regression analysis (solid black lines) for 
each family of mutations represents the $-value, as indicated, and 
can be compared with the slope for all other S4 mutations (solid 
gray lines). ANCOVA was used to determine statistical signifi- 
cance; *, P < 0.05. 



equilibrium constant for inactivation at a single voltage, 
normally 0-mV (-Ke q ,o = knaa.o/ &rec,o) , unidirectional rate 
constants at each voltage were derived by extrapolation 
from the linear portion of the relevant arm of the chev- 
ron plot (Fig. 1 D, i, dashed lines). After the introduc- 
tion of a point mutation, V535S, the plot of log(i^bs) 
against voltage demonstrates a clear depolarizing shift 
in the voltage dependence of inactivation, which is indi- 
cated by a slowing of kma.v and a concomitant accelera- 
tion of &rec,v- The resulting logarithm of the equilibrium 
constant for inactivation is clearly shifted more negative 
with respect to WT (Alog(i^ q , 0 ) = -0.91 ± 0.1; n=1), 
allowing a $-value to be calculated according to Eq. 7 
(see Materials and methods) . The kinetic properties of 
Kvll.l channel inactivation can therefore be exploited 
to allow its study using REFER analysis. 

To obtain an accurate estimate of a <&-value from a 
REFER plot, one or more mutations of a given residue 
must cause a sufficient perturbation to the transition 
pathway. We, both previously (Wang et al., 2011) and in 
Fig. S2, and others (Cymes et al., 2002; Fersht and Sato, 
2004) have established a Alog(,rv q ,o) of greater than 
or equal to ±0.5 as a reasonable cutoff criterion. Based 
on these criteria, we used REFER analysis to perform 
an extensive investigation of the role of the S4 helix 
during Kvll.l channel inactivation. To optimize our 
chances of finding mutant channels within the S4 helix 
that were amenable to REFER analysis, we performed 
three separate mutagenesis scans on residues Gly522 to 
Lys538. Each residue was mutated to alanine (Fig. 2 A), 
serine (Fig. 2 B), or tryptophan (Fig. 2 C). From these 
scans, several residues stood out as having at least one 
mutation with a Alog(^ q o) of greater than ±0.5 log units, 
namely, L529A/S, L530S/W, R531A/S, L532S, V535A/ 
W/S, and R537W (Fig. 2, closed bars, and Table SI). 
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Figure 5. Combinations of S4 mutants do not have additive ef- 
fects on the energetics of Kvll.l channel inactivation. (A) Bar 
graph comparing shifts in log(A^ qj o), relative to WT, for triple 
S4 mutant (L529S + L532S + V535S), double mutants L532S + 
V535S and L529S + V535S, and single mutants V535S, L532S, 
and L529S. Data are presented as means + SEM for 5-18 cells. 
(B) REFER plot of triple S4 mutant (L529S + L532S + V535S) , double 
mutants L532S + V535S and L529S + V535S, and single mutants 
V535S, L532S, and L529S. Data are presented as means ± SEM for 
5-18 cells. The linear slope (solid black line) indicates a $-value 
of 0.42 and was not significantly different to that derived from all 
of the individual S4 mutants (<£ = 0.50; P > 0.05 using ANCOVA) . 
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Although it is possible to derive ^-values from the in- 
dividual mutant channels listed above (see Table SI), it 
is important to note that the temporal information they 
provide is based on the effect of the specific mutation in 
perturbing the transition pathway, and may not neces- 
sarily reflect the role of the native amino acid residue 
per se. For example, any single mutation could affect 
two or more different processes that, by coincidence, 
cancel each other out to give a ^-value between 0 and 1 
(see Materials and methods) . However, this type of co- 
incidental combination will not occur with every muta- 
tion at that position. Therefore, a better estimation of 
€>-values can be derived from plots of log(&i, lact ,o) versus 
log(-Ke q ,o) for a family of mutations at a single residue or 
region of the channel (Grosman et al., 2000) . We there- 
fore interrogated S4 helix residues identified from our 
scanning mutagenesis in more detail by obtaining fami- 
lies of mutations at each position (Fig. 3) . Mutation 
of charged residues Arg531 (to Gin, Thr, and Asn) or 
Arg537 (to Gin, Asp, and Thr) did not produce any fur- 
ther mutant channels with adequate perturbations to 
Alog(i^ q>0 ) (Wang etal., 2011, and Fig. 3 A). The relative 
tolerance of neutral or polar side chains at these posi- 
tions, as well as at Lys525, Arg528, Arg534, and Lys538 
(Table SI and Wang et al., 2011), suggests that the posi- 
tively charged side chains are not critical to the inacti- 
vation gating transition. We therefore concentrated 
on the hydrophobic residues: Leu529, Leu530, Leu532, 
and Val535, which were further mutated to threonine, 
histidine, proline, asparagine, and glutamine. In the 
case of Val535, we also mutated this residue to glycine, 
leucine, isoleucine, methionine, and tyrosine. Perturba- 
tions to Alog(-Ke q ,o) of greater than ±0.5 log units were 
observed for several mutant channels at all positions, 
with the exception of Leu532 (Fig. 3 B) . Hydrophobic 
residues Leu530 and Val535 were particularly sensitive 
to mutation; i.e., tyrosine, methionine, glycine, proline, 
histidine, asparagine, and threonine mutations at Val535 
all exhibited a Alog(i^e q o) of >0.5 log units compared 
with WT, whereas only the relatively conserved hydro- 
phobic side chains of leucine and isoleucine produced 
smaller perturbations (Fig. 3 B). The relative intoler- 
ance of these hydrophobic residues to mutations that 
do not maintain their hydrophobic side chain proper- 
ties indicates their importance in mediating open-to- 
inactivated transition. 

The equilibrium constant represents the energy dif- 
ference between the open state and the inactivated state 
(see Fig. SI). Because energy and structure are related, 
a mutation-induced perturbation to the energetics of 
gating, observed as a shift in log(i^. q o), indicates that 
the amino acid residue experiences a change in envi- 
ronment at a given point during the transition from the 
open state to the inactivated state. The ^-value repre- 
sents the relative time point during the transition at 
which this change in environment occurs, with <l>-values 



close to 1 representing a change in environment early 
in the transition pathway, whereas a lvalue close to 0 
represents a change in environment at a time point 
closer to the end of the transition. To determine the 
relative timing of motion for each of the S4 helix hydro- 
phobic residues during the inactivation gating transi- 
tion, ^-values were calculated from families of mutations 
at each position. Plots of log(&i mc t,o) versus log(i^ q o) for 
all mutants at Leu529 (lvalue = 0.48 ±0.10), Leu530 
(lvalue = 0.33 + 0.07), Leu532 (lvalue = 0.68 + 0.08), 
and Val535 (€>-value = 0.49 ± 0.02) are shown in Fig. 4 A 
(i-iv) . In addition, a plot of all Leu529, Leu530, Leu532, 
and Val535 mutants is shown in Fig. 4 B (<&-value = 0.48 ± 
0.03) , and a plot of all mutants in the S4 helix is shown 
in Fig. 4 C (lvalue = 0.50 ± 0.02). Of the individual 
hydrophobic residues, only Leu530 produced a <I>-value 
that was significandy different from the <fr-value derived 
from the remainder of the S4 mutations (0.33 ± 0.07 for 
Leu530 vs. 0.52 + 0.02 for all other S4 mutants; analy- 
sis of covariance [ANCOVA]; P < 0.05), suggesting that 
Leu530 may be involved at a slightly later stage in the 
open-to-inactivated transition pathway. 

The observation that the slopes of the REFER plots 
for hydrophobic residues Leu529, Leu532, and Val535 
are not significantly different (^0.5) suggests that this 
region of the S4 helix moves as a single entity. To test 
this, we investigated whether combining V535S with ei- 
ther L529S and/ or L532S would give additive (expected 
if motions were independent) or nonadditive (expected 
for rigid body motions) effects on Alog^^o) values. 
Fig. 5 A shows a bar graph comparing the shifts in 
log( J K eq , 0 ) of single mutants (V535S, L532S, and L529S) , 
double mutants (L529S + V535S and L532S + V535S), 



r30 




Figure 6. Perturbations to inactivation gating and activation gat- 
ing show poor correlation for S4 residue mutations. Plot com- 
paring perturbations to AG 0 , versus WT, of inactivation gating 
(AAG° Inact , represented by Alog(_K,, q0 )) compared with perturba- 
tions to AG 0 , versus WT, of activation gating (AAG°A c ti V ; see Materi- 
als and methods) for each S4 residue mutation. Solid black line 
represents linear regression analysis constrained to go through 
WT, whereas dashed gray line is an unconstrained fit (R 2 = 0.38). 
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and the triple mutant (L529S+L532S+V535S; Triple 
S4) . Although the Alog(i^ qj0 ) was larger for the L529S + 
V535S double mutant when compared with the single 
mutations, the increase in Alog(^e q o) was not as large as 
the sum of the two individual residues combined. Fur- 
thermore, neither the double mutant L532S + V535S 
nor the triple S4 mutant produced Alog(.rv q ,o) values 
that were significantly different from the largest individual 
single mutant, V535S (Fig. 5 A). A plot of log(^ n act,o) 
versus log(.rv qj o) for all combination mutants produced 
a slope of 0.42 ± 0.07 and was not significantly different 
from that derived from the S4 helix as a whole (Fig. 5 B) . 
Our data strongly suggests that the face of the S4 helix 
that encompasses the hydrophobic residues Leu529, 
Leu532, and Val535 moves as a single entity during inac- 
tivation gating of Kvll.l channels. 
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Figure 7. Only combinations of S4 mutants alter the voltage 
sensitivity of Kvll.l channel inactivation. (A and B) Dot plots 
showing the linear slope values calculated from the derived uni- 
directional rate constants (shown as black lines in Fig. 1 D) for 
the onset of (hnna, A) and recovery from (k:^; B) inactivation for 
groups of mutants (to Ala, Ser, and Trp) at each residue in the 
S4 domain, as well as the double and triple mutants. Slope val- 
ues indicate the voltage dependence of the forward (inactivation; A) 
and reverse (recovery from inactivation; B) gating transitions. 
The mean of all mutations at each individual residue falls within 
±2 standard deviations (gray box) of the mean calculated from all 
single S4 mutants. 



Inactivation of Kvll.l channels is thought to be in- 
trinsically voltage dependent; i.e., the voltage depen- 
dence of inactivation is not directly linked to the voltage 
dependence of activation (Vandenberg et al., 2012). 
There are two ways that a mutation can affect the volt- 
age-dependent inactivation of Kvll.l channels. The first 



S4 vs S4-S5L 




S4 vs S5 



-10 12 
Inactivated [ J 

I Mi 




-10 12 

log (K 



eq,0' 




Open 



Figure 8. Multi-domain model of inactivation gating in Kvll.l 
channels. (A-D) REFER plots to compare S4 residue mutations 
(shown in black) against S4S5 linker (A), S5 (B), S5P linker (C), 
and S6 (D) residue mutations (shown in gray). ANCOVA shows 
that the linear regression slope, which represents the <l>-value, was 
significantly different for the S4 helix compared with the S5, S5P, 
and S6 helices (P < 0.05), but not the S4-S5 linker. (E) Cartoon 
showing two opposing subunits of Kvl 1 . 1 . Segments in the right- 
hand subunit are labeled S1-S6 and pore helix (PH). Segments 
in the left-hand subunit are color and number coded, according 
to derived <5-values (present study and Wang et al., 2011; Perry 
et al., 2013), to show the relative sequence of events that occur 
during the open-to-inactivated gating transition, where red (1) 
indicates the first step and purple (8) indicates the putative final 
step. The S4 helix (green, 5) , as well as the S4—S5 linker, expe- 
rience a change in environment shortly after the S5 helix and 
S5P linker but before the S6 helix. We propose that hydrophobic 
interactions between the S4 helix and the S5 helix help couple 
the voltage-sensor domain to the pore domain during the open- 
to-inactivated gating transition of Kvll.l channels. 
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is by altering the chemical free energy change associated 
with inactivation, i.e., a shift in the equilibrium set point 
of the voltage dependence (compared with WT), as re- 
flected by a change in log(,Keq,o)- The vast majority of 
point mutations in the S4 helix (66 out of 70; see Fig. 6 
and Table SI), and throughout the channel protein 
(Wang et al., 2011), caused a negative shift in log(i^ q o), 
reflecting a depolarizing shift in the voltage dependence 
of inactivation (i.e., stabilization of the open state). In 
contrast, the same set of S4 mutations altered the chem- 
ical free energy, and hence the voltage dependence, of 
activation gating in either direction to an approximately 
equal degree (Fig. 6). Thus, the S4 mutation-induced 
perturbations to the voltage dependence of inactivation 
show only a weak correlation (R 2 = 0.38) with the pertur- 
bations to the voltage dependence of activation (Fig. 6) , 
consistent with the notion that the activation and inacti- 
vation gating processes are not directly linked. 

The second way that a mutation can perturb inac- 
tivation is by directly altering the degree of voltage de- 
pendence, which can be detected as a change in the 
slope of the unidirectional forward (inactivation) or re- 
verse (recovery from inactivation) transition rate con- 
stants (shown as black lines in Fig. 1 D). Slope values 
representing the voltage dependence of inactivation 
(Fig. 7 A) and recovery from inactivation (Fig. 7 B) for 
all individual S4 residue mutations (to Ala, Ser, and 
Trp) tested during scanning mutagenesis of the S4 do- 
main all fall within ±2 standard deviations (gray box) of 
the mean. In other words, none of the individual residue 
mutations significantly alter the voltage dependence of 
inactivation gating. It was only when we combined hy- 
drophobic residue mutations (L532S + V535S for the 
forward rate) and (L529S + V535S and L529S + L532S + 
V535S for the reverse rate) that we observed modest 
changes to the voltage dependence of rates of inac- 
tivation (Fig. 7). Thus, it is possible that the S4 helix 
contributes to the voltage dependence of inactivation 
gating. However, as we only saw significant deviations 
with combined mutants, we should be cautious in sug- 
gesting that S4 plays a direct role. Even with this caveat, 



it is clear that the role of the S4 helix in inactivation gat- 
ing is different than the role it plays as the principal 
voltage sensor for activation gating in Kvll.l (Subbiah 
et al., 2004; Zhang et al., 2004; Piper et al., 2005) and 
other Kv channels (Swartz, 2008). 

REFER analysis demonstrates that hydrophobic resi- 
dues within the S4 helix are important determinants for 
inactivation gating in Kvl 1.1 channels. The relative tim- 
ing of the change in environment experienced by the 

54 hydrophobic residues during the inactivation gating 
transition can be determined by comparing the overall 
•Jj-value of M3.5 for the S4 helix with <J>-values derived 
from other domains within the channel protein (Wang 
et al., 2011), as shown in Fig. 8. The S4 <fr-value was not 
significantly different from that derived from the inter- 
nal S4-S5 linker (3> = 0.43 ± 0.04; P > 0.05; ANCOVA), 
suggesting that these two domains are coupled during 
inactivation gating. However, the change in environ- 
ment experienced by the S4 hydrophobic residues oc- 
curs just after a conformational change in the S5 helix 
($ = 0.78 ± 0.07; P < 0.05) and S5P linker (<t> = 0.60 + 
0.03; P < 0.05) but before a conformational change in 
the S6 helix (3> = 0.33 ± 0.05; P < 0.05) (Fig. 8) . 

When viewed on a Kvll.l channel homology model, 
based on the crystal structure of a Kvl .2/2.1 channel 
chimera (Long et al., 2007), the side chains of the hy- 
drophobic residues Leu529, Leu532, and Val535 face 
toward hydrophobic residues on the S5 helix of the 
neighboring subunit (Fig. 9) . This raises the question of 
whether the S4 helix and S5 helix are energetically 
coupled via an intersubunit hydrophobic interaction. 
To experimentally test for an energetic coupling between 
hydrophobic residues on the S4 helix (Leu529, Leu530, 
Leu532, or Val535) and hydrophobic residues on the 

55 helix (Ile560, Leu564, or Ile567), or a nonhydropho- 
bic residue on the S5P linker (Asp591), we combined 
perturbing mutations of residues within each helix 
(Fig. 10). If an energetic coupling exists between two 
residues, we would predict that the perturbations caused 
by the individual mutations, measured as Alog(i£e qj o) 
relative to WT, would not be additive when combined 
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Figure 9. Homology model of the Kvll.l channel. View 
parallel to the membrane of an amplified region of a Kvl 1 . 1 
channel homology model created using the Kvl. 2/2.1 chi- 
mera crystal structure (Long et al., 2007) as a template, ac- 
cording to the alignment shown in Fig. S3. The amplified 
region is indicated by the boxed region of the entire four- 
subunit homology model shown in the inset. Hydrophobic 
residues on the S4 helix (Leu529, Leu532, and Val535) of 
one subunit (subl; shown in green, with residues colored 
orange) face toward hydrophobic residues on the S5 helix 
(Ile560, Leu564, and Ile567) of the neighboring subunit 
(sub2; shown in blue, with residues colored purple). The 
arrow represents a cavity filled by lipid in the Kvl. 2/2.1 
crystal structure. 
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Figure 1 0. Double mutant cycle analysis of S4 and S5 helix resi- 
dues. (A) Schematic showing the principle of double mutant 
cycle analysis to test for an energetic coupling between Leu529 
on the S4 helix and Ile560 on the S5 helix. Perturbations to inac- 
tivation, measured by changes in log(X^ q> o) relative to WT, caused 
by individual mutations, L529S (AAG mml ) and I560A (AAG mut2 ), 
are compared with the perturbation caused by the double mutant 
L529S + I560A (AAG mu ,i +mut2 ) . Measured Alog(_Ke q?0 ) values repre- 
senting AAG mut i, AAG mut2 , and AAG mutl+mut2 are shown in paren- 
theses. An energetic coupling between Leu529 (S4) and Ile560 
(S5) residues is indicated by AAG mu ,i +mut2 being significantly dif- 
ferent from AAG mut i + AAG mut2 (P < 0.05; see Materials and meth- 
ods). (B-E) Double mutant cycle analysis for S4 helix mutations 
L529S (B), L530S (C), L532S (D), and V535S (E), paired with 



in the double mutant (i.e., AAG mut i +mut 2 < AAG mu ti + 
AAG mut 2; see Materials and methods) . Conversely, if the 
perturbations caused by the two single mutants were ad- 
ditive when combined in the double mutant, the two 
residues are not energetically coupled (i.e., AAG mut i +m ut2 = 
AAGmuti + AAG mut 2). For example, the schematic in 
Fig. 10 A represents double mutant cycle analysis be- 
tween two residues: Leu529 on the S4 helix and Ile560 
on the S5 helix. It is clear that the perturbations caused 
by the two mutations, L529S (AAG mutl ) and I560A 
(AAG mut 2), are not additive when combined in the dou- 
ble mutant, L529S + I560A (AAG mut i +mu t 2 ). In fact, the 
AAG nm ti+mut2 of the double mutant is significantly less 
than the theoretical additive value (AAG mut i + AAG mut2 ) 
(Fig. 10 B; P < 0.05; two-tailed Welch's t test) , indicating 
that Leu529 on the S4 helix and Ile560 on the S5 helix 
are, at least in part, energetically coupled. Similarly, com- 
bining the same S4 helix mutation, L529S, with a differ- 
ent S5 helix mutation, I567A, resulted in a AAG mut i +mut 2, 
which was significantly less than AAG mut i + AAG mut 9 
(Fig. 10 B; P < 0.05; two-tailed Welch's t test), suggesting 
that these two residues are also energetically coupled. 
In contrast, combining L529S with either L564A (S5 helix) 
or D591K (S5P linker) resulted in a AAG mut i +mut 2 ~ 
AAG mut i + AAG mut2 (Fig. 10 B; P = NS; two-tailed Welch's 
t test) , suggesting that Leu529 is not energetically cou- 
pled with either Leu564 (S5) orAsp591 (S5P). 

In addition to L529S, we also tested three other S4 
helix hydrophobic residue mutations: L530S (Fig. 10 C), 
L532S (Fig. 10 D), and V535S (Fig. 10 E). All show a 
similar pattern, in that the combinations with the S5 
residue mutations I560A or I567A are nonadditive (P < 
0.05; two-tailed Welch's t test), whereas the combina- 
tions with the S5 residue mutation L564A, or S5P linker 
residue mutation D591K, are additive. One exception 
to this is the combination of L532S (S4) with D591K 
(S5P linker), which resulted in a nonadditive effect. Al- 
though we cannot rule out the possibility that these two 
residues are energetically coupled, it seems unlikely 
given their distant spatial proximity in models of the 
Kvll.l channel (Tseng et al., 2007; Durdagi etal., 2012). 
An alternative explanation is that the combination of 
mutations, L532S + D591K, introduces a non-native in- 
teraction or alters the structure of either or both of the 
two stable end states. Although this is not indicated by an 
abnormal lvalue ($ = 0.75 + 0.01; n = 10; see Table SI) , 
that is to say outside of the range 0 and 1, it remains 
possible that the double mutant (L532S + D591K) affects 
two or more processes that coincidentally cancel each 
other out to give a cfr-value between 0 and 1. 



mutations on the S5 helix (I560A, L564A, or I567A) or S5P linker 
(D591K). Measured values of Alog(Ae q0 ) for individual or double 
mutants are shown as closed bars, whereas the theoretical addi- 
tive values, AAG mutl + AAG mut2 , are shown by open bars. Data pre- 
sented as means ± SEM for 6-14 cells (see Table SI). 
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In general, however, the nonadditive effects of com- 
bining S4 helix mutations (L529S, L530S, L532S, or 
V535S) with S5 helix mutations (I560A or I567A) pro- 
vide strong evidence for the notion of hydrophobic 
coupling between the S4 helix and the S5 helix during 
inactivation gating in Kvll.l channels. Based on the 
Kvll.l channel homology model, we suggest that this 
energetic coupling is likely to be caused by an intersub- 
unit interaction, which is critical for inactivation gating 
in Kvll.l channels. 

DISCUSSION 

C-type inactivation gating in potassium channels is thought 
to culminate in a conformational rearrangement of the 
selectivity filter, which results in the loss of one or more 
K + coordination sites (Cuello et al., 2010). Recendy, we 
used REFER analysis to show that inactivation of Kvll.l 
channels is a complex phenomenon in which a pre- 
cisely timed sequence of conformational changes in 
multiple interconnected domains of the channel must 
occur before reaching the nonconducting stable end 
state, a process analogous to the opening of a Japanese 
puzzle box (Fig. 8) (Wang et al., 2011). In the present 
study, we used REFER analysis to examine the role of the 
S4 "voltage-sensor" helix during inactivation gating. Our 
results indicate that in contrast to channel opening, where 
the positively charged residues within the S4 helix are cru- 
cial for "sensing" changes in membrane potential to trig- 
ger opening of the S6 activation gate (Subbiah et al., 2004; 
Zhang et al., 2004; Piper et al., 2005), it is predominantly 
the hydrophobic residues within the S4 helix that are the 
crucial molecular determinants for inactivation gating. An 
overall lvalue of M).5 indicates that the S4 helix, together 
with the internal S4— S5 linker, experiences a change in 
environment after conformational changes in the S5 helix 
(3> = 0.78) and S5P linker (€> = 0.60) but before a confor- 
mational change in the S6 helix (€> = 0.3) (Fig. 8). Based 
on our Kvll.l homology model and double mutant cycle 
analysis data, we propose that intersubunit interactions 
between hydrophobic residues within the S4 and S5 heli- 
ces are important determinants of inactivation gating in 
Kvll.l channels. 

Although we observed perturbations to the energetics 
of inactivation gating upon mutation of both charged 
and hydrophobic residues within the S4 helix (Fig. 2) , it 
is clear that the charged residues exhibit greater toler- 
ance to nonconserved side chain replacements than do 
the hydrophobic residues, at least at positions 530 and 
535 (Fig. 3) . This finding is initially surprising given the 
crucial role of the charged residues in mediating the 
voltage dependence of Kvll.l channel opening (Subbiah 
et al., 2004; Zhang et al., 2004; Piper et al., 2005). The 
voltage-sensor domain, however, is likely to adopt at 
least three stable conformations: resting, active, and re- 
laxed (Villalba-Galea et al., 2008) . The transition from 



the resting to active conformations is linked to opening 
of the intracellular channel gate (Swartz, 2008). Our 
voltage protocols were designed to ensure that the S4 
was predominantly in the "active" position before mea- 
suring the onset and recovery from inactivation. In this 
context, it is possible that the hydrophobic S4 residues 
may contribute to so-called voltage-sensor relaxation. 
Whether this in turn is coupled to inactivation we can- 
not determine from our studies and is a source of con- 
troversy (Villalba-Galea et al., 2008; Lacroix et al., 201 1; 
Shirokov, 2011). Regardless of whether voltage-sensor 
relaxation is involved or not, what is clear from our ex- 
periments is that the hydrophobic S4 residues are key 
determinants of the open-to-inactivated transition in 
Kvll.l channels. 

The voltage dependence of activation and inactiva- 
tion gating in Kvll.l channels are clearly separable 
(Wang et al., 1997; Zou etal., 1998; Johnson et al., 1999; 
Sanguinetti and Xu, 1999; Piper et al., 2005; Clarke 
et al., 2006; Vandenberg et al., 2006). Our data indicate a 
poor correlation between the mutation-induced changes 
in the equilibrium of inactivation gating compared with 
that of activation gating (Fig. 6) , supporting the notion 
that the S4 involvement during these two gating transi- 
tions involves a quite distinct molecular mechanism. 
Furthermore, none of the individual S4 residues we 
mutated direcdy perturbed the voltage dependence of 
either the forward or the reverse inactivation gating tran- 
sitions (Fig. 7) , suggesting that the origin of the voltage 
dependence of inactivation may lie beyond the region 
of the S4 helix. 

REFER analysis provides information regarding the 
relative time point at which a perturbation affects the 
native transition pathway of a reaction (Auerbach, 
2007) and can be used to infer when the mutated resi- 
due or region of the protein was involved in that transi- 
tion pathway. Previously (Wang et al. , 20 1 1 ) , we observed 
that altering the external concentration of potassium 
ions gave a <I , -value of ^1, suggesting that a loss of po- 
tassium ions from the selectivity filter is the first step 
during the inactivation gating transition. Diminishing 
•^-values were then observed for families of mutations 
in the S5 helix (€> = 0.75), S5P linker (3> = 0.6), S4-S5 
linker ($ = 0.45), S6 helix (3> = 0.3), and pore helix 
(<!> = 0.25). In this context, we can confirm that the S4 
helix (<1> = 0.50) experiences a change in environment 
after the S5 helix and S5P linker but before the S6 helix 
(Fig. 8) . An overall lvalue of ^0.5 for the S4 helix was 
similar to that derived from the internal S4-S5 linker, 
suggesting that these two domains are coupled during 
inactivation gating. 

Combining three of the S4 hydrophobic residue mu- 
tations (Leu529, Leu532, and Val535) produced shifts 
in log(Xe q o) that were clearly not additive, suggesting 
that these residues experience a change in environ- 
ment at a similar time point. Families of mutations at 
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Leu529, Leu532, or Val535 gave <l>-values that were not 
significantly different from that derived from the re- 
mainder of the S4 domain, whereas the ^-value derived 
from Leu530 mutations was significantly smaller (€> = 
M).3) . This suggests that Leu530 may undergo a change 
in energetics at a time point closer to that of the S6 do- 
main (<I> = 0.3) than the other S4 residues. According to 
our homology model (Fig. 9), Leu530 is not in direct 
contact with the S6 helix but instead faces the S3 helix 
and/or the lipid membrane. It appears unlikely, there- 
fore, that Leu530 directly interacts with S6 helix residues, 
but it is entirely possible that the two regions of the chan- 
nel undergo independent conformational changes at 
approximately the same time point during the transi- 
tion pathway. Nevertheless, we suggest that the relatively 
consistent nature of the ^-values that we derived from 
mutations of residues throughout the S4 helix indicates 
that the majority of the S4 domain, along with the inter- 
nal S4— S5 helix, experiences a change in environment at 
the same relative time point during the open-to-inactivated 
gating transition. 

Although REFER analysis provides temporal informa- 
tion about the role of the S4 helix during inactivation 
gating, it cannot provide a direct indication of the na- 
ture of the change in environment experienced by the 
S4 helix. We therefore used a combination of homology 
modeling and double mutant cycle analysis to probe 
the nature of this structural change. Based on our 
Kvll.l channel homology model (Fig. 9), we hypothe- 
sized that the important hydrophobic residues in the 
S4 helix face directly toward the S5 helix of the neigh- 
boring subunit. Previously, several different types of in- 
teractions between the voltage-sensor domain and the 
pore domain have been suggested in various subtypes 
of potassium channel. In Kv channels, including Kvl 1.1, 
residues at the intracellular end of the S4 helix, to- 
gether with the S4-S5 linker, interact with residues at 
the intracellular end of the S6 helices that form the ac- 
tivation gate (Lu et al., 2002; Tristani-Firouzi et al., 
2002; Long et al, 2005; Ferrer et al, 2006). This type of 
interaction is important in coupling the transition of 
the voltage sensor into the active conformation upon 
membrane depolarization to the opening of the internal 
S6 activation gate (Long et al., 2005). A second set of 
interactions has been proposed between residues at the 
extracellular ends of the SI (Long et al., 2005, 2007; Lee 
et al., 2009; Bocksteins et al., 2011) or S4 helices (Elinder 
et al., 2001; Neale et al., 2003; Bocksteins et al., 2011) 
with residues at the extracellular end of the S5 helix of 
the neighboring subunit. It has been suggested that 
these interactions could help constrain the voltage sen- 
sor to the pore during channel opening (Lee et al., 2009) . 
Lastly, in hyperpolarization-activated RATI channels, 
there is evidence for interactions between two non- 
charged S4 residues and two residues in the middle of 
the S5 helix (Lai et al., 2005). 



Our double mutant cycle analysis data support the no- 
tion of an energetic coupling between hydrophobic resi- 
dues on the S4 helix (Leu529, Leu530, Leu532, and 
Val535) and hydrophobic residues on the S5 helix (Ile560 
and Ile567) of the neighboring subunit during inactiva- 
tion gating. However, the number of different side chain 
mutations tolerated at some of these positions, namely 
529 and 532, suggests that these two regions are not 
tightly packed. Leu529 through Val535 lie within the in- 
tracellular half of the S4 helix and, as such, are below the 
region proposed for the extracellular voltage sensor- 
pore domain interaction. They are also clearly distinct to 
the residues at the intracellular end of S4 that have been 
suggested to form interactions with the intracellular end 
of S6, i.e., Asp540 with Leu666 (Ferrer et al., 2006) . They 
could, however, be consistent with the midlevel interac- 
tion observed in RATI channels (Lai et al., 2005). Based 
on our homology model, it is possible that some of the 
hydrophobic residues are close enough for a direct hy- 
drophobic interaction, such as Leu529 (S4) with Ile567 
(S5), but in other instances, the distance between resi- 
dues does not support such a notion. In the Rvl. 2/2.1 
chimera x-ray structure, there is a gap between this re- 
gion of the S4 helix and the S5 helix of the neighboring 
subunit that is occupied by lipid (Long et al., 2007). It is 
conceivable that the replacement of bulky hydrophobic 
residues with serine affects the interactions between the 
hydrophobic residues in S4 and S5 with the hydrophobic 
lipid. The S4 and S5 domains may then be energetically 
coupled through a shared lipid interaction. It has been 
shown that gating of Rv channels is influenced by the 
lipid composition of membranes (Schmidt et al., 2006, 
2009) . In the case of RvAP channels, Schmidt et al. (2009) 
showed that the lipid composition influenced multiple 
gating transitions, with the most dramatic effect observed 
for rates of recovery from inactivation. Investigation of 
this hypothesis in Rvl 1.1 channels will, however, be diffi- 
cult, as it would require expression of sufficiendy large 
amounts of recombinant protein to allow reconstitution 
into different lipid bilayers. 

Most voltage-gated channels contain several pairs of 
hydrophobic residues, interspersed between positively 
charged residues, along the length of the S4 segment 
(Anselmi et al., 2007; Long et al., 2007; Pless et al., 
2011). The highly conserved nature of the S4 hydro- 
phobic residues across all Rv channels (Anselmi et al., 
2007; Long et al., 2007; Pless et al., 2011) suggests that 
the S4 and S5 hydrophobic interaction we propose 
could be common to all members of the family, rather 
than being unique to Rvl 1.1 channels. In fact, it has 
been suggested that small changes to the degree of hy- 
drophobicity within the S4 helix could alter the extent 
of packing against the S5 helix, and that this could re- 
flect differences in the kinetics of voltage-sensor motion 
during channel opening (Anselmi et al., 2007). How- 
ever, it should be noted that the sequence homology 
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and hydrophobicity of the S5 helix is much less con- 
served across the Kv family (Shealy et al., 2003; Ju et al., 
2009) , and that C-type inactivation gating in other Kv 
channels is not intrinsically voltage dependent (Hoshi 
and Armstrong, 2013). Because inactivation gating has 
much slower kinetics than those of channel opening in 
other Kv channels, it may be difficult to separate these 
two gating processes sufficiendy to test whether hydro- 
phobic residues within the S4 domain play a similar role 
in mediating selectivity filter gating in other Kv chan- 
nels. Nevertheless, it is clear that S4 hydrophobic resi- 
dues are critical mediators of inactivation gating in 
Kvll.l channels. 
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